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I. INTRODUCTION 



Since quark spin accounts for only a small fraction, E « 0.3, of the nucleon's spin H'a one 
wonders where the rest of the spin resides. Sehgal II pointed out long ago that experimental 
data on hyperon /3-decays and the assumption that any strange quarks in the nucleon are 
unpolarized o requires £ ~ 0.6 and he suggested that quark orbital angular momentum 
{Lq) - expected in relativistic quark models - was a likely candidate. 

More recently, the possibility that gluons may carry a significant fraction of the nucleon's 
spin has been raised and debated, ua Much debate has centered on the possibility that a 
gluonic contribution may "contaminate" the axial charge sum rules used to extract £ from 
polarized lepton scattering data, leading some to speculate that T is large and positive. 
Whatever their impact on axial charges, there is no doubt that gluon spin (T) and orbital 
angular momentum {Lq) can contribute to the nucleon spin. In QCD the nucleon's total 
angular momentum can be written as 

l - = L Q + L G + h: + Y. (i) 

The definition of £ and its relation to the triangle anomaly were explored shortly after the 
publication of the SMC data. @~& However, the definitions of the other terms and their 
relation to local operators have received less attention. No one knows how to measure Lq 
and L G . Y is defined as the integral of the polarized gluon distribution, 

T{Q 2 ) = [ 1 dxAg{x,Q 2 ) = f 1 dx{g^{x,Q 2 )-g l {x,Q 2 )) (2) 
Jo Jo 

in the sense of a parton model sum rule. 

The purpose of this Letter is to study T. First we review the definition of Ag in terms of 
gauge invariant operators. Next we integrate Ag, and specialize to A + = gauge, where 
T can be related to the matrix element of a local gluon operator. We discuss the gauge 
dependence of T and its relation to the spin generators obtained via Noether's theorem. We 
conclude that T is given by the matrix element of products of gluon vector potentials and 
field strengths in the nucleon rest frame and in A + = gauge, 

T(Q 2 ) = ^(e 3 | 2 Tr{{E x A) 3 + A ± ■ B±}\ Q2 \h) (3) 

where _L refers to the directions transverse to the e% axis defined by the target spin, and Q 2 
is the renormalization scale of the operators.^ 

Eq. (H) is amenable to evaluation in constituent quark models. Just as quark models 
make predictions for E which are tied to the phenomenology of /3-decay axial currents, so 
too they make predictions for T, tied to the phenomenology of baryon mass differences. 
Quark models claim that baryon spin splittings {e.g. — M^) originate in lowest order 



An earlier version of this paper contained an incorrect version of this equation which was incor- 
rectly claimed to be gauge invariant. The author wishes to thank Ian Balitsky and Xiangdong Ji 
for raising questions of gauge invariance and for helpful discussions on other issues. 
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exchange of transverse (magnetic) gluons. The predictions agree well with experiment, and 
the presentations have found their way into textbooks. Known as "color magnetism" , 
this effect implicitly requires ambient color magnetic fields within hadrons. The model 
calculations are performed to lowest non-trivial order in olqcd, where the gluons behave like 
eight Abelian vector fields coupled to color. The spin dependent piece of the Born graph for 
gluon exchange between quarks i and j can be rewritten as 

8 r 

AM = ~ E E / d 3 xB?(x) ■ B°(x) (4) 

i>j a=l 

just as in electrodynamics. Thus the baryons contain color magnetic (and electric) fields 
which may carry spin angular momentum. The sign of T is correlated with the sign of baryon 
spin splittings and originates in the non-abelian character of QCD forces. If quark spin 
forces were abelian, the A would be lighter than the nucleon and T would be positive. The 
magnitude of T depends on details of quark model wavefunctions and on the renormaliation 
scale assigned to the quark model calculations. However the model's prediction of the A — N 
mass difference constrains the size of T to some extent. 

A generic non-relativistic quark model (NQM) yields Tnqm = ~ f °^f M (^)i the standard 
bag model yields r& ag = — 0.1a;, ag . Note that 2T is the gluon spin fraction of the nucleon 
spin. We hesitate to assign a specific value to a for fear the resulting number will be 
taken too seriously. On the other hand the model parameters are constrained to some 
extent by the magnitude of baryon spin splittings making a numerical estimate possible. 
In the non-relativistic quark model the N — A splitting is AM = 8n ^2 M (<5 3 (r))- If ; fo r 
example, we choose Gaussian wave-functions scaled to reproduce the proton's charge radius, 
and a quark mass m q ps Mjy/3 to obtain approximately correct magnetic moments, then 
AM ps 0.35q;jvqmA^v, whence «atqm ~ 0.9. With this value of cxnqm we find Ttvqm ~ —0.8. 
In the bag model, cibag ~ 2 in order to fit baryon mass differences, whence r& aff ps —0.2. 

Note that these estimates apply at the quark model renormalization scale. To obtain a 
prediction for Q 2 ~ 2 — 10 GeV 2 relevant to experiment, it is necessary consider the evolution 
of T with renormalization scale. It has long been known that T evolves homogeneously, 
and that the quantity a(Q 2 )T(Q 2 ) is a renormalization group invariant to leading order. O 
Therefore the sign of T should be a reliable prediction of the models. If we take the numerical 
estimates seriously, evolution will tend to bring the two predictions closer together because 
the bag estimate will be multiplied by a larger factor of a{^)/a{Q 2 ). For example, at 
ctQCD ~ 1, ^bag ~ -0.4, and T NQM ps -0.7. 

In the next section we introduce and analyze the operator measure of gluon spin. We 
discuss gauge invariance. We show that the operator we evaluate in the rest frame is the 
same one measured by integrating the gluon helicity asymmetry Ag(x,Q 2 ). In Section 3, 
we estimate the magnitude of T in simple models. In Section 4 we discuss the reliability of 
our calculation and and mention further applications. 



II. THE OPERATOR DESCRIPTION OF GLUON SPIN 

Experimenters will not measure ((E x A) 3 + A± ■ B±) directly. Instead they will measure 
the polarized gluon distribution function, Ag(x, Q 2 ), in deep inelastic lepton scattering. Ag 
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measures the probability to find a gluon with its helicity parallel to the nucleoli's helicity 
minus the probability to find it antiparallel. The integral of Ag measures V and therefore 
will tell us which local operator, if any, is to be associated with the gluon spin. We begin 
with the operator representation of the polarized gluon distribution function, H 

Ag(x,Q 2 ) = J rfre- MrP+ (P,e 3 |Tr{F +a (r,0)J(r,0)F a + (0)}| Q2 |P,e 3 ) 

+ (x — > —x), (5) 

where F (and A) are matrices (F = X)o=i -^ a A a , etc.). {F a } are in the adjoint, and {A a } 
are in the triplet representation (with Tr{A a A 6 } = \5 ab ). ^ ± ,^ _L are light-cone coordinates 
and (£~, 0) denotes the point £~, £ + = = 0. The label Q 2 is a reminder that the tower of 
local matrix elements in the Taylor expansion of FF are understood to be renormalized at 
a factorization scale, Q 2 , and finally X is the Wilson-line integral, 

J(r , 0) = Pexp Ug J* dy-A + (y-, 0)) . (6) 

The standard parton interpretation follows from eq. (^) if we choose A + = gauge and 
introduce the momentum decomposition of the fields F and F quantized at £ + = 0. 

In order to integrate eq. (^) over x we must study the apparent singularity at x — 0. 
Physically, Ag(a;,(5 2 ) is not expected to diverge as fast as 1/x, so the ^"-integral must 
vanish as x — > 0. This means that the singularity at x = is integrable. If Ag were found 
to diverge like 1/x or faster, our analysis would have to be reconsidered. With this in mind, 
we can interchange the x and £~ integrations and, because Ag is symmetric in x, we can 
use the principal value prescription at x = 0, -^^£ e - lcra = _ iite(ot). We obtain 

T(Q 2 ) = 2^+/ ^r£(r)(^e 3 |Tr{F+«(r)X(r,0)F Q +(0)}| Q2 |P,e 3 }. (7) 

This expression cannot be simplified further unless we choose A + = gauge. 

In A + = gauge X = 1 and P +a = ^|^v4 a , so we may perform the £~ integration. The 
terms at £~ = ±oo vanish because the integral in eq. (§]) converges when x — 0. Only the 
local (£~ = 0) contribution survives. We choose the rest frame for P, and after some algebra 
we are left with the expression we seek, 

T(Q 2 ) = -^(P, e 3 |2 Tr {^P+ 2 - A 2 F +1 }| q JP, e 3 > 

1 ■ (P, e 3 |2 Tr {(E X I) 3 + 1 ± • P ± } IP, e 3 >, (8) 
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where P l = P*°, and P* = —^e^ k F^ k . The choice of A + = gauge was essential to this 
derivation. Without it, T does not appear to be associated with a local operator. 

A + = does not completely fix the gauge in QCD. Residual non-abelian gauge transfor- 
mations, SA^ = d^5a + [5a, A 11 ], are allowed provided they obey d + 5a = which preserves 
A + = 0. T must be invariant under this residual gauge symmetry. Using the Bianchi identity 
- [D*, F uX ] + [D U , F x »} + [P>\ F^} = — it is straightforward to show that 
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5T oc (P, e 3 |Tr {ad + F 12 }\P, e 3 > (9) 

which vanishes because a is independent of £~ and the + derivative sits inside a forward 
matrix element. Thus T is invariant under the residual gauge symmetries of A + = gauge. 

Eq. (H) is not a familiar representation for the spin angular momentum stored in a 
gauge vector field. Happily it can be related directly to the gluon spin term in the angular 
momentum tensor density in QCD. Angular momentum in QCD, as in any field theory, is 
described by a rank-3 Lorentz tensor, M fluX . In a physical gauge, where there are no ghosts, 

MSod = \^ (x x d» - x»d x ) i/, + h^ lal5 ip 

- 2 Tr {F^ a (x u d x - x x d u ) A a } + 2 Tr {F» X A V + F V »A X } 



- - Tr F 2 (afgi* - x x g' xv ) , (10) 

The second term in eq. ([10"|) measures the quark spin - at least up to subtleties arising from 
the triangle anomaly. HiS The first and third terms look like the quark and gluon orbital 
angular momentum respectively, because they have the standard "connective" form of orbital 
angular momentum in a field theory, W(x x V)$, where II and $ are canonical coordinate 
and momentum respectively. The last term contributes only to boosts. The fourth term is 
a candidate for the gluon spin. Let us define, 

Mf x = 2 Tr {F flv A x + F^A V }. (11) 

Then comparison with eq. (|8]) shows that 

nQ 2 ) = 7^-(P,e 3 \M+ 12 \P,e 3 ), (12) 
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in A + = gauge. This identification makes physical sense since the parton model distribu- 
tion should measure helicity along the e 3 -axis (hence v — 1, A = 2) in an infinite momentum 
frame, which corresponds to \i — + in the laboratory. The restriction to A + = gauge 
is natural in the parton model. Needless to say, this restriction does not render F gauge 
dependent: There is a corresponding operator definition of T in any gauge. However it is 
not simple or even local. Note that this definition of gluon spin does not correspond to the 
gluon piece of the generator of rotations in the laboratory, which would be M° 12 and would 
not naturally appear in light-cone gauge. 

This discussion suggests that a natural definition of the quark and gluon orbital angular 
momentum might select the +12 component of the appropriate piece of M >luX in eq. 
While this is attractive, a physically interesting Lq or Lq must be determined by what 
can be measured experimentally. As long as no measure of orbital angular momentum is 
experimentally accessible, the definition will remain open to question. 



2 Eq. (|io|) was derived in Ref. Q, where the generator of gluon spin rotations, was incorrectly 
identified as A x E. Since this sign is crucial to our results we have checked a) that canonical 
tranformations give E x A and b) that the two gluon contributions combine to give a total angular 
momentum ofJ = fx£x8as expected. 
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III. THE GLUON SPIN IN QUARK MODELS 



Quark models of the light hadrons fall into two general classes: non-relativistic quark 
models, where quarks are described by the Schroedinger equation (perhaps including rela- 
tivistic corrections) and confined by some two body color dependent forces; and bag models, 
where relativistic quarks, governed by the Dirac equation, move in some confining back- 
ground field imagined to be self-consistently generated by their deformation of the non- 
perturbative QCD vacuum. Both extremes give good explanations of the mass spectrum of 
the lightest hadrons (pseudoscalar and vector mesons, and octet and decuplet baryons). A 
major role is played by color mediated, spin dependent forces. 00 In this Section we shall 
see that the gluons responsible for these spin splittings are anti-aligned with the nucleon spin 
(r < 0); that this is a particular consequence of the non-abelian nature of QCD interactions, 
and that the effect has roughly the same magnitude in both types of model. 

The quantity we want to evaluate in models is[] 

r(^) = (T,e 3 \ J d 3 x2Tr{(E(x) x A(x)) 3 + A ± (x) ■ B ± (x)}\T, e 3 > (13) 

One way to evaluate this expression would be to compute the relevant Feynman diagram 
- obtained by inserting this gluonic operator into the lowest order Born diagram for gluon 
exchange between bound quarks (in A + = gauge). The graphical method would require 
us to construct and use confined gluon Green's functions in A + = gauge, which is unnec- 
essarily complicated. It is easier to compute the ambient color fields by directly integrating 
the QCD equations of motion (which reduce to eight copies of Maxwell's equations in the 
Abelian approximation). In this way we obtain expressions for A, B, and E which depend 
on the quark color, spin and spatial coordinates and can be evaluated with the help of model 
wavef unctions. 

In a generic quark model the color-electric fields will be given by the gradient of a time- 
independent function of the quark degrees of freedom, E a (x) = — V$ a (x), with 

3 a (*) = |:EW,2<) (14) 

i 

$ is an operator in the space of the quark color (A") and position (xi) states. It depends 
on some model Green's function Q . For example, in an unconfined, non-relativistic model 
Q = l/\x-Xi\. 

The magnetic field, B a (x) can likewise be written as the curl of a time independent 
function of the quark variables, B a (x) = V x U a (x), with 

U a (^) = fj2^iX0(x,x i ). (15) 



Quark model states are more conveniently normalized to unity than covariantly, so we change 
normalization accordingly, (P\O(0)\P) = 2M J d 3 x{T\0(x)\T) , where |T) is a quark model state 
normalized to unity. 
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For example, Q = (x — Xi)/2m q \x — Xj| 3 in a non-relativistic model (where m" = g<Ji\% /2m q 
is the quark's color magnetic moment operator). The nucleon polarization (and the gauge 
choice, A + = 0) selects the 623-axis. Matrix elements of the operator di will therefore vanish 
except in the e 3 direction. So we conclude that ([/" ■ e 3 ) = (£/f 3 ) = 0. 

So far this result is quite general. It holds in any quark model where the gluons are 
treated to lowest order (abelian approximation) and U has no component along the nucleon 
spin. In "symmetric" quark models the quark distributions in the nucleon ground state 
are not correlated with the overall spin. When integrated over the wavefunction of the j th 
quark, then, the resulting color electric field is radial and independent of j, 

- J cPxjVg&ZjMixj)] 2 = ^Q(r), (16) 

where r = \x\ and Q(r) = An /J* dr'r' 2 \ip(r')\ 2 is the color charge inside the sphere with radius 
r. Likewise integration over xi simplifies Q, 

J (PxiG&x^Xi)^ = xh(r), (17) 

where h(r) describes the vector potential generated by the model dependent magnetization 
density. 

The operators $ and U are not yet the appropriate scalar and vector potentials for 
the gluon field because they do not satisfy the A + = gauge condition. However suitable 
potentials are easily constructed. Define 

A 0a (x) = $ a (f) 

A a (x ) = U a (x ) - V f d($ a (x, y, C) (18) 
Jo 

These potentials generate E a and B a in the usual way and satisfy the gauge constraint, 
A 0a + A ia = (remembering U 3a = 0). The choice of time independent potentials as well 
as the lower limit on the ( integration correspond to residual gauge freedom available in 
A + = gauge. 

We now substitute the operator definitions of A 0a and A a into eq. ( |T3| ) and obtain 

+ U a ±l (x) ■ B a ±j (x) + V ± f?(x) ■ B a ±j (x)}\T,e 3 ) (19) 

where we have separated out the contributions from individual quarks % and j to each of 
the field operators. Note we have dropped the i = j terms. These correspond to "self- 
angular momentum" effects and are associated with questions of renormalization, which are 
discussed further in Section 4. 

The second and third terms in eq. flT9| ) vanish. The second vanishes after integration by 
parts because V x £° = 0. There is no associated surface term. The third term vanishes 
for spatially symmetric quark wavef unctions. 

Consider now the fourth term and write out the space components explicitly (suppressing 
the color (a) and quark (i, j) labels, and bras and kets, 



E?{x) x [//(f) + Et(x) x V//(f 
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^idfdU, dfdU 2 



\dy dz dx dz J 
= J d 3 x (E 1 U 2 -E 2 U 1 ) 

+ j s d 2 se 3 ■ f (u^x) £ d(E 2 (x, y, () - U 2 (x) jT d(E 1 (x, y, C)) (20) 

where the surface integral is over a sphere at large distance (for unconfined gluons) or the 
bag surface (for bag-like models). The first term in eq. (|20|) is identical to the first term in 
eq. 

We combine these results and substitute the parameterizations of E and U from eq's. ([16]) 
and QT7D to obtain, 

r = la I* drrQ{r){h{R) - 2h(r)), (21) 
9 Jo 

where the r-integration goes to infinity in generic non-relativistic quark models, but ends at 
r = R, the bag surface in the bag model. The term proportional to —2h(r) is the volume 
integral of E x U, and the h(R)-term is the surface contribution left over from integration 
by parts. In reaching eq. (pl|) we have used Z)o=i K = ~ 2/3 for i ^ j. [Because the 
nucleon is a color singlet, [J2a=i Sf=i ^j] 2 — 0. The result follows since Sa =1 (A") 2 = 4/3.] 
Also, a = g 2 /4:iT and J^i^ji^f) = 2 in a state polarized along the e 3 -direction. 

Now let us specialize to quark potential models where quarks are confined but color 
is not. Gluon field strengths fall off at large distances like abelian multipoles. The non- 
relativistic vector potential, U = fa x r/r 3 , corresponds to h(r) nqm oc l/2m q r 3 . At large 
R, Q(R) —> 1, so the surface term in eq. fl2ip vanishes. Substituting for Q(r) in terms of 
the quark wavefunction, if)(r) and interchanging integrations, eq. ( |2~TD reduces to 

8 r°° 
FnqaM) = —^at N Q M (i4) J o drr\i)(r)\ 2 

8 1 
= -^r- a NQM{^ )(-), (22) 

where we have restored the quark model renormalization scale, to remind us that this 
value pertains to some low scale at which the model is formulated. 

The parameters m q , (-) and Ozvqm^o) are an m °del dependent, but not unconstrained. 
m q pa 0.3GeV reproduces nucleon magnetic moments. Another constraint comes from the 
A — N mass difference which is given by 



AM =- aNQM ^ l) (5\r)) (23) 
3 m 1 



in the non-relativistic quark model. To obtain a numerical estimate, we assume a gaussian 
wavefunction adjusted to reproduce the root-mean-square charge radius of the proton. Then 
we find a NQM (fJ,l) pa 0.9, and T NQM (n%) pa -0.8. 

Static bag model calculations of hadron spinsplittings were carried out by explicit con- 
struction of color electric and magnetic fields. We may borrow results from that work 
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to evaluate eq. (|21). The color magnetic field is calculated from the QCD generalization of 
Maxwell's equation, V x B a = gip^a\ a ip augmented by the boundary condition f x B a = 
at r = R. A short calculation yields, 



pr 

Q(r) = / 
J o 

Hr) = {\ 

rr 

nif) = 

J o 




rfrV 2 (/ 2 (r') + g 2 (r')) 



dr'—r 3 f(r)g(r) 



[ dr' ^ ^ r ^ }, where 



(24) 



where /(r) oc jo(xor/R) and g(r) oc ji(xor/R). [xo is the lowest solution to the eigenvalue 
condition tana; = x/1 — x (xq = 2.0428).] Substituting explicit wavefunctions we find 



Standard bag model calculations of baryon spin splittings require olqcd ~ 2, so we find 
Tfeag ~ —0.2 at the renormalization scale of the model. 



Our calculations force us to conclude that the gluons responsible for the spin splittings 
among light baryons are anti-aligned with the spin of their parent nucleon. This exacerbates 
rather than helps resolve the problem with the nucleon spin, suggesting we look elsewhere for 
a large positive contribution to the spin. Our calculation required two assumptions beyond 
the traditional formulation of quark models: first, that "self-interaction" contributions to 
T could be ignored, and second, that the estimate of T at the quark model scale, /i 2 , has 
something to do with its value at experimentally accessible scales. Here we will first look at 
the reliability of our estimate in the context of the models. Then we will comment further 
on the two assumptions and finally mention some possible extensions of the work. 

The troubling sign of T does not appear to depend on the details of the models. The 
directions of the color electric and magnetic fields are fixed by the charges and currents 
which give successful descriptions of a variety of hadronic phenomena including magnetic 
moments as well as colormagnetic spin splittings. The minus sign arises because the gluons 
are non-abelian — the two spectators of any given quark in the nucleon generate color fields 
appropriate to an antiquark. The same effect flips the sign of baryon spin splittings. If 
quark spin forces were abelian, T would be positive, but the A would be lighter than the 
nucleon. 

On the other hand, the magnitude is of T quite uncertain. The value of T(Q 2 ) at scales 
relevant to experiment varies quadratically with the assumed value of a(/ig) taken from 
fits to baryon spin splittings. One factor of a comes from the operator itself; the other 
factor arises in scaling from ^ to Q 2 . There is reason to believe that symmmetric quark 
models overestimate the value of a necessary to account for spin splittings because they 
ignore the correlations color magnetic interactions would introduce into wavefunctions. We 
expect that the color magnetism would correlate quark pairs into color 3, spin diquarks 
(the most attractive channel) with the effect of increasing the matrix element and decreasing 




(25) 



IV. DISCUSSION AND CONCLUSIONS 
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the value of a necessary to reproduce spin splittings. An estimate of the size of this effect 
goes beyond the scope of this Letter. However we caution the reader against taking the 
calculated magnitude of T very seriously. The most important conclusion of our work is the 
prediction that the sign of T should be negative. If experiment finds T to be positive, and if 
neither self-interactions nor evolution change the model predictions, then the quark model 
calculations of baryon spin splittings will have to be re-evaluated. 

Regarding the self interaction terms, the proper prescription in principle is "calculate all 
graphs and renormalize" . Consider the case of an isolated electron in QED: the E and B 
fields which surround it contribute a log-divergent term to the expectation of Mjt 12 . The 
same occurs for the stress tensor (one piece of which is the self-energy). Renormalization 
removes these infinities. Since M^ u \ = x^yx — xJT^x (up to a total derivative) renormalizing 
T also renormalizes M. If we renormalize on-shell, the spin is found on the renormalized 
electron line and the self-fields are to be ignored. The same argument would seem to carry 
over to the quark model and dictate that we ignore all terms with % — j in eq. (0). However 
we cannot renormalize a confined quark on mass-shell, so we cannot assume that all self- 
field effects can be renormalized away. The problem also arises in the bag model where the 
calculation and renormalization of the self-field contribution to Mjt 12 would yield a finite 
result analogous to the Lamb shift correction to the self-energy of a bound electron. This 
program has actually been carried out for the self-energy of a quark in a bag. @ Only a 
complete treatment of the renormalization of a confined quark would enable one to decide 
how much of the self-field should be included in the normalization of a quark at the scale 
/z 2 , and how much should be attributed to its binding inside the nucleon. We have, in effect, 
assumed that all self-field effects are to be included in the definition of the renormalized 
quark operators and omitted from our calculation of T. This assumption should be examined 
more closely; it is possible that finite parts of the i = j terms make a significant contribution 
to T at the scale /i 2 ,. 

The renormalization point and scheme dependence of T has been discussed extensively. 
In A + = gauge it is easy to see that the forward matrix element of the Kogut-Susskind 
current, coincides with our definition of T. So the studies of the QCD evolution 

of the forward matrix elements of (K + ) IE apply without alteration to the evolution of 
T. To leading order aQCD(Q 2 )F(Q 2 ) is independent of Q 2 . Beyond leading order the Q 2 
dependence of T is scheme dependent, and not of much interest until one understands the 
renormalization scale and scheme dependence of quark model calculations. These problems 
notwithstanding, it appears that the quark model result T < is preserved by leading order 
evolution. 

The analysis presented here raises several questions and suggests some avenues for future 
study. First, no doubt, is to understand better the origin of the result T^qm < and 
its generality. The light-cone perturbation theory formalism pioneered by Brodsky and 
Lepage, 111 and applied to the parton structure of positronium by Burkardt, @ seems like 
the natural framework in which to explore the result further. Another direction would be to 
use these models to explore the x-dependence of the nucleon's polarized gluon distribution 
(at the scale /i 2 ,), since the shape of Ag is of interest. These calculations might serve to 
guide the application of sophisticated (next-to-leading order) perturbative QCD fits to spin 
dependent gluon distributions. §1 Finally we note that there is much less reason to compute 
the spin- independent gluon distribution using these models. The lowest moment of g(x, Q 2 ) 
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does not converge and is not associated with a local operator. Furthermore models that take 
seriously the gluon role in confinement, leave us puzzled whether the long-range confining 
field (e.g. the bag itself in the bag model) should be associated with a (spin-independent) 
gluon distribution or not. 
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